Localisation of fermions to brane: Codimension d > 2 

W. Nahm*, and D. H. Tchrakian^* 

*School of Theoretical Physics, Dubhn Institute for Advanced Studies 
10 Burhngton Road, Dubhn 4, Ireland 

^Department of Mathematical Physics, National University of Ireland Maynooth, 

Maynooth, Ireland 

February 1, 2008 



Abstract 

We investigate 4 + d dimensional fermionic models in which the system in codimension-d 
supports a topologically stable solution, and in which the fermion may be localised to the brane, 
with power law in 'instanton' backgrounds and exponentially in 'soliton' backgrounds. When 
the fermions are isoscalars, the mechanism fails, while for isospinor fermions it is successful. 
As backgrounds we consider instantons of Yang-Mills and sigma models in even codimensions, 
solitons of sigma models in odd codimensions, as well as solitons of Higgs and Goldstone models 
in all codimensions. 



1 Introduction 



The idea that our world is a (fairly flat) 3-brane in a higher dimensional space has deep roots in the 
19th century |2H]- It was introduced into physics in the context of cosmological defects on one hand 
and the branes of M-theory on the other. In any case one has to explain why the observed fermion 
masses are so much smaller than the mass scale given by the transversal extension of the brane. 
On some parameter range one expects the branes to approach classical configurations, such that a 
semiclassical description of matter on the brane is possible. Qualitative features like the appearance 
of low energy excitations on the brane should be independent of the parameters at least locally, such 
that these semiclassical descriptions may be relevant for realistic cases. 

We study the extension of the original work of Rubakov and Shaposhnikov PP localising a fermion 
to the brane, in which a 5 dimensional model, i.e. one with codimension-1 is employed. There ^1 
121, the brane Lagrangian is the (p'^ (or sine-Gordon) system in 1 dimension, supporting the kink- 
soliton. We are concerned with the possibilty of extending this mechanism to the case of arbitrary 
codimension-c? > 2. In general, fermions are localised when their reduced Dirac equation in the d 
dimensions transversal to the brane has zero modes. This always is the case when the index of the 
corresponding Dirac operator is positive. 

In case [1^ the wave function of the fermion drops exponentially away from the brane. Such 
backgrounds will be called solitonic. For them an explicit formula for the index is known [3], and 
physically relevant specialisations are known both in odd and in even dimension |31 13] . The index 
only depends on the behaviour of the scalar fields at infinity, gauge fields are assumed to approach 
zero at infinity and do not influence the index in a direct way. Nevertheless, the corresponding 
internal symmetries are crucial. In the backgrounds of Yang-Mills or sigma model instantons, the 
fermionic fields have a power-law fall-off away from the brane. In such cases the available index 
calculations depend on compactifications at infinity, which allow an application of the Atiyah-Singer 
index theorem 0. 

For solitonic backgrounds, non-zero indices occur in both even and odd dimensions, whereas in 
the instanton case one needs even dimensions. Even when the index is zero or negative, solutions of 
the reduced Dirac equation might exist. For fermions in instanton backgrounds this can be excluded 
by vanishing theorems for the kernel of the corresponding Dirac operator (the Weizenbock formula), 
but in other contexts little is known. In particular, no general results for solitonic backgrounds seem 
to exist. 

The index formulas require cumbersome if straightforward integrations, but often an analytic 
approach is possible. When a background can be deformed into a superposition of well separated 
configurations with spherical symmetry, the index is additive and only these spherically symmetric 
situations have to be considered. For any of them the index is the sum over the indices for fixed 
angular momentum, which are determined by the behaviour of ordinary differential equations and 
easy to calculate. We consider some families of such configurations. As new features we obtain that 
certain solutions are universal and apply to any number of transversal dimensions, and in certain 
cases a vanishing theorem can be established. 

We have considered two alternative types of fermionic models in 4 + o? dimensions. In the first 
of these, (i), the Dirac spinor is isoscalar, such that the Dirac operator on the codimension involves 
only unmodified partial derivatives and no gauge field. In the second type, (ii), the Dirac spinor on 
the codimension is in general an isospinor under an internal SO{d) group, i.e. it is a square matrix 
valued array. The Dirac operator features a covariant derivative, but the index is determined for the 
case of a vanishing gauge field. Here we also consider an exceptional d = 2 case when the covariant 
derivative is Abelian and the Dirac spinor is a two-component column. 
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Like in [H 12] , our brane Lagrangians in d dimensions support finite energy topologically stable 
solutions, both of 'instanton' and of 'soliton' types. Our nomenclature here is the following: Consider 
the topological charge density g[(p], e.g. the Chern-Pontryagin (C-P) density and its descendents 
when [ip] symbolises YM fields (in even dimensions) and YM-Higgs fields, respectively, or, the degree 
of the map when [ip] stands for sigma model fields. The scalar deansity g[p] is essentially total 
divergence, in the sense that its variation vanishes when it is subjected to arbitrary variations 6ip. 
(Indeed in the (C-P) case g[ip] is precisely a total divergence Q[(p] = diQi, Qi being the Chern-Simons 
density.) When g[ip] is subjected to spherical symmetry, which is the case for the fields [p] at infinity, 
it reduces always to a total derivative of some function a, i.e. 



We denote this function pertaining to the classical solutions pc by (Xc = (^['^c]- The profiles (fT^ 
and (fT3jl stated below in section 2.1 specify our nomenclature of 'instanton' and 'soliton' types, 
respectively. Thus, finite energy topolologically stable solutions to Yang-Mills (YM) or a sigma 
model systems in even d are typified by 'instanton' profiles, (|Tl|l . and solutions to sigma models in 
odd d, or to Higgs models, i.e. YM-Higgs (YMH) systems, or Goldstone models, in both even and 
odd d, by 'soliton' profiles (fTK|l . By Goldstone model, in turn, we mean the gauge decoupled version 
of a Higgs model, provided of course that the soliton persists in the gauge decoupling limit of the 
YMH model in question. In this context, the solutions of the symmetry breaking model(s) employed 
in P3I2] are typified by 'soliton', (fT3)). profiles. 

The main feature of the mechanism of jUEj is dimensional descent of a 4+(i dimensional fermionic 
model to 4 dimensional Minkowski space, the Dirac field on which is assumed to be chiral and to be 
massless. For codimension-ci > 2, after the descent one is left with a nontrivial Dirac equation on the 
codimension-d, which we have called the residual Dirac equation. The crucial step is that of finding 
normalisable zero modes of this residual Dirac equation. The asymptotic behaviour of the residual 
Dirac spinor is then responsible for the decay of the Dirac field off the brane. Whether this decay is 
achieved, and if so, is it power like or exponential, will depend on the A + d dimensional fermionic 
model chosen. 

It will turn out that the desired normalisable zero modes do not exist for models of type (z) for 
any codimension d, and that such solutions exist for models of type (ii) for all d. 

To date, extension to higher codimensions {d > 2) for this mechanism has been performed for 
codimension-2, in a series of works started by Libanov and Troitsky [3 IE] in flat space, and in the 
presence of gravity, in p. The models employed in both j2j and jHI are of the type {ii), namely 
featuring a covariant derivative in the Dirac operator. In these models [3 IH] , what guarantees the 
existence of zero modes of the residual Dirac equation in codimension-2, is the choice of such a Yukawa 
coupling in the 6 dimensional model that leads to a residual Dirac equation which coincides with the 
particular Dirac equation on an Abelian-Higgs background for which Jackiw and Rossi ^Ul have 
constructed the zero modes explicitly. To follow this line of approach for codimensions-c? > 3, one is 
naturally led to considering such 4 + d dimensional fermionic models which result in d dimensional 
residual Dirac equations for which we know there exist normalisable zero modes, or better still that we 
can construct such solutions. Dirac equations in d dimensions in the background of a YM 'instanton' 
or a YMH 'soliton' (or its associated Goldstone 'soliton') supporting such zero modes are the natural 
candidates which will be proposed. Dirac equations in d dimensions will be solved in the appropriate 
'instanton' or 'soliton' backgrounds, supported respectively by a hierarchy of YM and YMH models 
in these dimensions. 

The first in the hierarchy of YM models is the usual YM system in 4 dimensions, and its extensions 




(1) 
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to all even dimensions as given in ^T] support 'instantons' in these dimensions, analogous to the 
4 dimensional instanton 17. It is these 'instantons' in d = 2n dimensions that we will employ as 
backgrounds for the construction of the zero modes of the residual Dirac equations, extending^ the 
d = 4 result of Jackiw and Rebbi jil3j to arbitrary d. 

The first in the hierarchy of non Abelian YMH models is the Georgi-Glashow (Higgs) model in 
d = 3, which supports the the usual monopole |15]. The zero modes of the Dirac equation on this 
background was given long ago by Jackiw and Rebbi ^Hl- The extension to arbitrary d is completely 
straightforward: The result of ^H] follows directly from the existence of monopoles |il5j in the 3 
dimensional Higgs model. The corresponding solitons of Higgs models in arbitrary dimensions d have 
been systematically shown to exist, being constructed numerically in jT71^1|T^. The generalisation 
of the d = 3 result of jTHI to arbitary d then follows almost trivially. Moreover, the adaptation of the 
result of to the case of Goldstone model backgrounds also follows systematically, by employing 
the solitons presented in [201 IHl 122] • 

In section 2 we present the i + d dimensional models on the space with coordinates xm = {x^, Xm), 
/i = 0, .., 3 labeling the Minkowski space and m = 1, ..2,d the codimension, and, we give the Ansatz 
separating the varaibles and Xm in the field equations. This describes the dimensional descent. 
Both type {i) and type (ii) models, i.e. with Dirac operators featuring both partial derivatives 
and covariant derivatives, are presented in this section. The resulting residual Dirac equations in d 
dimensions will then be examined in detail in the subsequent sections 3 and 4 respectively. In 3, 
it will be shown that for models of type (i) the fermion cannot be localised to the brane for any 
d. In section 4 type (ii) models will be analysed. In the first subsection of 4, zero modes of the 
residual Dirac equations in even d > 4 dimensional 'instanton' backgrounds of YM systems will be 
constructed, resulting in the power localisation of the fermion to the brane. In the second subsection 
of 4, the corresponding zero modes in all ci > 2 dimensional 'soliton' backgrounds of Higgs (or their 
associated Goldstone) systems will be constructed, resulting in the exponential localisation of the 
fermion to the brane. A summary of the results is given in section 5, and three appendices have been 
supplied. Appendix A describes the (even) d = 2n dimensional YM models and their 'instantons'. 
Appendix B describes the d dimensional Higgs models and their 'solitons', and Appendix C describes 
the Goldstone counterparts of the latter. 

2 The model(s) and residual Dirac equations 

We will consider the following two types of fermionic actions, formally expressed as 



The first of these, (0), pertains to the type (i) family of models featuring partial derivatives in the 
codimension. The components of the spinor field on the codimension-rf, in are isoscalar. a{ip\ 
in (j2)) is a scalar function of the fields [if] symbolising the scalar and/or the YM field describing the 
brane Lagrangian in the codimension-ci. Specifically, it will be defined as the leading term in the 
spherically symmetric restriction of the topological current, e.g. the Chern-Simons term in the case 
of YM. Both 'instanton' and 'soliton' backgrounds can be accommodated in this scheme. 

^We restrict our considerations to Dirac equations on spherically symmetric backgrounds only and exclude the 
multicentre backgrounds employed in or even periodic backgrounds used in |14|. 




(3) 



(2) 
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The second, represents type (ii) models in which Dm = {d^,Dm)- The components of the 
spinor field on the codimension-d, in (jH)), are isospinor except in the case d = 2 case when they 
are isoscalar. These models subdivide further into two subclasses, ones with /i = 0, i.e. without a 
Yukawa term, and, those with 7^ 0. The models Q with /i = accommodate only 'instanton' 
backgrounds, while those with /i 7^ accommodate only 'soliton' backgrounds, and in this case 
is a matrix valued function of [ip\. 

The Dirac equations arising from (j2)) and Q are, respectively 

V^dM^ + li(T[ip,]i! = (4) 

f^DM^ + /iH[(^,]^ = 0, (5) 

to be solved on the classical background ^d^^) to be precised later. In the present work, we anticipate 
the use of radially symmetric background solutions in terms of the radial variable r = \xm\ of the 
co-dimension, though a richer spectrum of such backgrounds arises when this symmetry is relaxed. 

Denoting the 4-dimensional (spacetime) coordinates by and the coordinates of the codimension 
by Xm, we represent the 4 + dimensional gamma matrices T m = (f ^, Tm) by 

f = 7m ® I > f m = 75 ® (6) 

in terms of the 4-dimensional gamma matrices 7^ and their chiral matrix 75, and the ci- dimensional 
gamma matrices r^. 

Our separability Ansatz, which also effects the dimensional descent, is 

^(X^, X^) = ^(X^) (g) ^{Xm) , (7) 

and when applicable, 

E[ip{Xm)] = ^®^HXm)]. (8) 

In (jHI), the array ^[(/^(xm)] is a matrix valued array whose size will be determined by the representaion 
in which the gauge connection in the covariant derivative Dm is. (In the generic case this will be 
the Gamma matrix representation of SO{d).) Detailed definitions of ^ for particular models to be 
considered, are given in section 2.2. 

Using (jni) and the separability Ansatz (III),®, the Dirac equations (jlj) and © yield, respectively 

7'^^^^ ® + 7^* ® ^"dm^P + /i a, ^ ® ^ = (9) 
7^9^* ® + 7^^ ® r™D„V^ + /X ^ ® ^ = 0, (10) 

in which we have used the notation cTc = cr[y9c] and = ^[vj- If we now invoke the existence of the 
zero modes of the Dirac field in 4 dimensional spacetime 

Yd^^ = , 

and require that the Dirac spinor is chiral, i.e. that 

= ^ , 

then (Pj). ()10|) finally reduce to the residual Dirac equations in d dimensions 

(ra^ + z^ajV = (11) 
{T^'Dm + I^Q^P = 0. (12) 
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It is in order to mention at this point, that for the case of /i 7^ models 0, the covariant derivative 
Dm in (fT^ will eventually be replaced by the partial derivative dm, since the detailed analysis of 
(fT^ to be carried out subsequently is restricted only to Goldstone models associated to the Higgs 
models described in Appendix C, namely to the gauge decoupled versions of the associated Higgs 
models described in the Appendix B. We stress that all our results are valid also for the Higgs 
model 'soliton' backgrounds, and the only reason we eschew working with the latter is that the 
corresponding analysis of (jl2p yields qualitatively the same results as in the Goldstone case, and the 
analysis in the latter case is somewhat simpler. 
We will seek solutions to (fTT |) .(fT ^ satisfying 

V''(0) < 00, ||V'||<oo. (13) 

It will turn out that there exist no solutions of (fTTj) satisfying (fT^ for any d, but will find such 
solutions to (fT^ for all d. 



2.1 Definitions of a[(p] 

a[(p] are defined as the leading terms in the topological currents at infinity, i.e. when the fields defining 
them are spherically symmetric. This subsection is subdivided in two parts. 

In the first, we give the definition of cr[(p] for a topologically stable background supported by a 
0{d + 1) sigma model, as well as a background supported by a SO{d = 2n) YM system in even, 
d = 2n, dimensions. The latter are the spherical-symmetrically restricted Chern-Simons densities 
of SO{d = 2n) YM field, which turn out to be described by essentially the same 'instanton' type 
profile 

-l^a, ^1, (14) 

r^O r— >oo 

as in the case of even dimensional sigma models. 

In the second we consider backgrounds supported by Higgs or Goldstone models, typified by 
'soliton' profiles 

^ a, ^ 1 , (15) 

r<— r— ►oo 

as in the case of odd dimensional sigma models. Detailed expositions of (jl4j) and (|15p are given in 
the following two subsections. 



2.1.1 a[ip\ for 0{d+ 1) sigma model and SO{d = 2n) pure YM backgrounds 

d dimensional 0{d + 1) sigma models and their topologically stable solitons have been discussed 
extensively elsewhere |22j so we do not elaborate on them here. Best known amongst these is the 
d = 2 dimenesional scale invariant O {3) sigma model whose solitons, namely the well known Belavin- 
Polyakov vortices pi] , are evaluated in closed form. Here we are concerned only with the topological 
boundary conditions the relevant solitons satisfy. Moreover, as noted above, we will restrict to the 
case of radial (spherically symmetric) solitons. So we state these, in terms of the d + 1 component 
scalar fields = ix^^X'^), subject to the constraint = 1: 

X"^ = sin/(r)x" , x'+'=cos/(r) (16) 

In ()16|) x"^ = r~^x^ is the unit vector in the codimension. The topological charges stabilising 
the solitons of these models are the winding numbers, which take on unit values provided that the 
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solutions satisfy the asymptotic conditions^ 

lim/(r) = , lim /(r) = vr , (17) 

r—*0 r^oo 

Our definition of the function = <^[f{r)] corresponding to the solution = fi^) of this model, is 
that given by ((T)). 

We list the functions <Jc{d) for this model, for the cases d = 2, d = 4 and d = 3, 5, separately for 
even and odd d. Up to unimportant multiplicative constant depending on the angular volumes, we 
find following P^ . 

a{d = 2) (X cos/ , a{d = 4) (X (cos/ — -cos^/ 
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for even d = 2,4 respectively, and 



a(rf = 3) oc /- ^sin2/ , a(ci = 5) oc ( |/ - sin 2/ + ^ sin4/ ) , ... (19) 



for odd d = 3,5 respectively. We see that the ranges of these topological charge densities are quite 
different for even and odd d. 

(fT8|l and (fT9|l result, qualitatively, in the following profiles of Uc 

+ 1< — ac — ^ -1 , for even d (20) 

r-s— r^oo 

0< — ac — ^ 1 , for odd ci . (21) 

r'<— 1 — >cxD 

We next turn to the SO{d) YM system in even (i-diemensions, for which the spherically symmetric 
Ansatz, analogous to (fT6|) . is 

_ l-w{r) _ _ 1 / i±rvi\ 

Analogously to ()17|1 . the 'instanton' boundary conditions that result in topologically stable (anti)- 
selfdual solutions to the systems of YM hierarchies [11], are 

limw(r) = +1 , lim w{r) = —1 , (23) 

r— >0 r— >oo 

which coincides with (pT|l under the replacement 

cos/(r)< ^ w{r) , ... (24) 

Now in all even diemensions there exist Chern-Pontryagin charge densities, whose spherically sym- 
metric restrictions are the analogues of (|18p. Up to unimportant numerical factors, these densities 
in dimensions d = 4 and d = 6 are 

a{d = 4) oc ^w- ^ly^^ , a{d = 6) oc - '^w^ + ^w^^ . (25) 

Note that the first {d = 4) member of (j25|l coincides with the second [d = 4) member of (|T8|) under 
the replacement ()24|1 . This is a recurring coincidence. It is obvious now that the profile of Uc in this 
case coincides with fEIH) . 

^The more usual alternative limr^o fir) = tt , \imr~^oo f{r) = is not adopted, for the sake of making contact 
with the usual asymptotics for the corresponding YM fields. 
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2.1.2 a[Lp\ for SO{d) Higgs/Goldstone model backgrounds 

We now consider topologically stable backgrounds supported by d dimensional SO{d) Higgs mod- 
els [HI [H [H] and their associated Goldstone models [201 12111221 • 

The topologiocal charges of these models are described by scalar fields (j)"^ , m = 1,2, ...,d, in d 
dimensions. Apart from the various kinetic terms, Goldstone models are distinguished by a symmetry 
breaking self-interaction potential, leading to the important asymptotic condition 

lim|0"^P = r72 (26) 

r— >oo 

in which 1] is the VEV with inverse dimensions of length. Here again, we restrict to the radially 
symmetric fields 

(f)"' = 7] h{r) x"" , (27) 
and for the special case oi d = 2 dimensions, the radially symmetric vorticity n field is 

(j)-^ = r]h{r)rf^ , n*" = (cosn(/., sinn</.) . (28) 

The topological charges stabilising the solitons of these models are the winding numbers, which take 
on the unit value for the following asymptotic conditions 

lim/i(r) = , lim h{r) = 1 . (29) 

r— >0 r^oo 

The function cTc is now expressed in terms of the classical soliton profile h{r). One difference from 
the sigma models of the previous subsection however is, that the d = 1 case for Goldstone models, 
unlike the sigma models, does not trivialise but coincides, for example, with the ip^ model. Another 
difference is that the winding number density for the radial fields (12 7|) does not take qualitatively 
different expressions for even and odd d, as in (fT5|) - (fT^ . 

In both even and odd ci-dimensional Higgs and Goldstone models, the leading term in the winding 
number density turns out to be proportional to 

h{ry a, =%-'^|0'"|'^ = //"V (30) 

and since in the following we will need only the asymptotic values and not detailed behaviours of 
(Jc(r), we omit the d-th power of h{r) in (jHUj) and simply state the topologically meaningful asymptotic 
behaviour 

0^^ae™l , for all d. (31) 

2.2 Definitions of 

Unlike the quantities a[(p] presented in the previous subsection, which are isoscalar, the quantities 
^[if] in (jHJ are matrices with isotopic indices. In models employing sigma model or YM 'instanton' 
backgrounds in even codimension-rf, /i = so that S[<y9] is defined only for models employing 'soliton' 
backgrounds with yU 7^ 0. With odd d sigma model backgrounds in turn, there is no useful efinition 
for ^[v?]. The reason is simple, and hinges on the requirement that the residual Dirac equation ()12j). 
like ()11|) . should develop a mass term asymptotically in the codimension. 

Let us examine the Yukawa term in ()12|) in the asymptotic region, which is subject to spherical 
symmetry. Consider the matrix valued function ^[ip] in (jH)) in terms of the two alternative codimension 
fields x"" = (x^^X*^^^)) pertaining to the 0{d+ 1) sigma model, and 0"^, to the Higgs or Goldstone 
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model. The only natural forms for matrix valued ^ in the gamma matrix represenation in isospace 
are proportional to the following matrix valued quantities 

ecxr^x" , eocr^0"^, (32) 

respectively. Inspecting the spherically symmetric Ansatze f|TH|l . P7j) and the asymptotics (fTTj) . 
()29|) required, one sees that the only Yukawa term which leads to a nonvanishing mass term is the 
second member of ()32|) . Henceforth, models of type (jH)) with yU > will be restricted to Higgs or 
Goldstone model 'soliton' type backgrounds only, with the corresponding Yukawa term determined 
by the second member of ((221) • The type (ii) models considered are typified by the definitions of the 
quantity ^. There will be two such choices. 

The first applies in the case where the gauge connection in the covariant derivative is Abelian, 
which is the case only for d = 2, e.g. in the background of the usual Abelian Higgs model. In this 
case, it is possible to take the Dirac field to be an isoscalar, and our choice for ^ is 

e = aicr^0™ = 0^11 + ^0^3 ,m = l,2, (33) 

which is (Euclidean) Lorentz invariant. 

The second concerns the case of generic codimension-ci, where the SO{d) connection is non 
Abelian, our choice for ^ is 

e = i®r„0"^, (34) 

where the matrix 1 is labeled by the spinor indices and the matrix Tm4>"^ is labeled by the isospinor 
indices. While fl34j) is defined for non Abelian backgrounds, i.e. for c? > 3, it applies also to the 
d = 2 case formally. In that case, we express the Abelian gauge connection, say am, in formally 
antihemitian form 

Am = - am CTs (35) 

acting on the matrix valued Higgs field $ = (f)"^am- The covariant derivative in (jl2j) for the d = 2 
case of the generic model is then defined by the connection 



3 Type (i) models with isoscalar 'ilj{m) 

We will show that type (z) models with action (0), on backgrounds with either type of profile (jl4p 
and (fT^ of the function dc, do not support solutions satisfying condition (fT^. for any codimension 
d. This section is divided into four subsections. In 3.1, 3.2 and 3.3 we analyse the residual Dirac 
equation (fTTj) for d = 2, d = 3 and arbitrary d respectively. This yields a set of coupled ordinary 
differential equations, which are of the same form for all d. Then in 3.4 we show that these equations 
do not have solutions satisfying the required property (fT^ . 

3.1 Codimension d = 2 

The 2 component residual spinor ip{xm) is subjeted to radial symmetry 

^=\fl e-'<^ ) ' (36) 
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with m and m', both integers. Denoting m,m' instead by for uniformity of notation for all d, 
the variables r and in equation (lllj) separate for I' = / + 1, resulting in the pair of coupled first 
order equations 

f[--fi + ^cf2 = (37a) 
r 

/2 + ^/2 + ae/l = 0. (37b) 

r 

3.2 Codimension d = 3 

The residual 2 component spinor ip{xm) transforms as a spin-| spinor under 3 dimensional rotations, 
and to achieve a separation of variables we employ the spinor harmonics [211 expand 

^ = h{r)n^^ + /2(r)4;;. (38) 
The spinor harmonics are defined as 

^t} = c{i, 1,1 ± If"' 1 y,, 1 (^, 0) x^ + c{i, 1,1 ± if-"'^ - 1 y,, {e, <p) x-i m 

in which Yi^rn{0, <P) are the spherical harmonics and x±i are the constant valued 2 component eigen- 
vectors for spin- 1 

X.J = (J) . X-i = {\). (40) 
Evaluating the Clebsch-Gordan coefficients in (jH^ and substituting (^01) . we have 



'+'"+^y,._i(^,0) \ -J^Y, 



The result of acting with the gradient operator on the spherical harmonics can be systematically 
calculated applying the Clebsch-Gordan series [26j. Applying this to the residual Dirac equation 
(fTTjl with the Ansatz pHjl . and setting /' = / + 1, we have 



(/( - ^ /i + f2^ nl^l^ + (/2 + ^ /2 + /i) nl^ = 



(42) 



resulting in 



- -fi = -CTJ2 (43a) 
r 

/2 + ^/2 = -crc/i. (43b) 
r 



The similarity of ()43a|) - ()43b|) with ()37a|) - ()37b)) is manifest and holds for arbitrary codimension, as we 
shall see immediately below. 
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3.3 Arbitrary codimension d 

To generalise to arbitrary codimension d, we need some preparation. Let pk = —idk, (Jij = i[Ti, Tj]/2, 
i, j = 1, . . . d, and 



1 " 

A. = — ^ ] aij(xiPj — XjPi). 



2 



Squaring yields = (ra — 2) A + L^, where 



XjPif. 



i<j 

It is well known that the eigenvalues of are l{l + d — 2), I = 0,1, ... . The corresponding possible 
eigenvalues of A are — / and l + d — 2, but by inspection only for / = 0. Thus the list of A eigenvalues 
is 0, —1, —2, . . . plus d — l,d,d+ 1 . . .. 

In analogy to the case = 3 we introduce the spinor harmonics f^l^f ' ^1m' where 

L'nlS = Kl + d-2)Q'^, 
•^^ Hm ~ " Hm ■ 

The index M stands collectively for eigenvalues of angular momentum operators which commute 
with A. Determining the spinor harmonics explicitly is more complicated than for d = 3, but not 
necessary for our purpose. 

The operator V^Xm transforms eigenspaces for — / and {d — 1 + I) into each other, since 

A{nj = {nj{d-i-A), 

as can be checked easily by multiplying the gamma matrices and using obvious symmetries. Thus 
we can put 

q(-) _ -pm - q(+) 

where {d — 2 + I') = [d — 1 + I), thus /' = Z + 1 as for d = 2,3. Conversely we have 
Another multiplication of gamma matrices yields the Dirac operator in the form 

Pm — r Xm ~~r~ H r a. 

dr r 

With 

the Dirac equation takes the form 

f[ - -fi + cjj2 = (44a) 
r 

f + l±±Zlf, + CTe/l = . (44b) 

r 

We see that both ipTj) and for d = 2 and d = 3 are of the same form as for arbitrary d. 
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3.4 Nonexistence 



The presentation here is adapted to both sigma model as well as to Goldstone model backgrounds, 
namely for both 'inatanton', (fT^ . and 'soliton', ()15|) . profiles of crdr). Let us analyse the d = 2 
equation ()37|) . noting that the same conclusions hold for the arbitrary case ()44|) . 

In the r -C 1 region, the solutions of equations ()37al37b|) which are differentiable at the origin, 
have the asymptotic forms 



In particular, /2(0) = and /i(0) is finite for all possible values of /. At infinity, both fi and /2 
decay exponentially. 

Now the equations (j37j) yield 



which clearly is impossible. There are therefore no solutions of ()12|1 satisfying ()18|l in both the 
backrounds (fT^ and (fTH|l . and in any d. 

4 Type (n) models 

Type (ii) models with action Q separate in two main subclasses, namely those with /i = 0, presented 
in the first subsection 4.1, and those with /i > 0, presented in the second subsection 4.2. The quantity 
E[(p{xm)], defined in section 2.2, will be specified further in two cases. The first of these is a particular 
model with codimension-c? = 2 leading to the recovery of the result of |71IHI, presented in subsection 
4.2.1, while the second is for the generic codimension-rf > 2 models, presented in subsection 4.2.2. 

/i = models of type (ii) are defined exclusively in even co dimension- c/(= 2n), since YM 'in- 
stanton' backgrounds exist only in even dimensions ^. fi > models of type (m) are defined for all 
codimension-(i, in Higgs or Goldstone 'soliton' backgrounds. In the latter case we shall eschew the 
Higgs backgrounds in favour of the corresponding associated Goldstone backgrounds in which the 
gauge field is suppressed. 

4.1 Type {ii) models with /i = 0: even > 4 isospinor ^(x^) 

Here, the solitons in whose background the residual Dirac equation ()12|) is to be solved are restricted 
to the 'instantons' of d = 2n dimensional YM models described in Appendix A. To solve equation 

^It is also possible for Grassmannian sigma models, in which case the Dirac operator in ^ features a composite 
connection in terms of the Grassmannian field. Once we know this is a possibility, it is superfluous to present it in 
detail here. 





-/l/{ + /2/^ + -[//l' + (/+l)/2'] =0. 



Integrating over r and using the said boundary conditions, one finds 
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(|12|) with /i = in the spherically symmetric background for the YM connection (j22j) . we subject 
the isospinor ipi^m) to spherical symmetry 

^ = /i(r)]l + /2(r)SWx„. (45) 

Substituting (j45|) and the radial Ansatz (|22j) for the YM connection into the residual Dirac equation 
()12|) . with = 0, we find the familiar solutions 



/2 = r-^('^-l)e-^('^-l)/^'^^ (46b) 

Adopting the asymptotics (j^ . it is easy to see that /i(r) satisfies the required condition (fT^ . while 
/2(r) does not and must be rejected^. It also follows from that the relevant solution, e.g. /i(r) 
has a power decay at infinity. 

4.2 Type (m) models with fi ^ 

In this case, the separability Ansatz (jH)), and the resulting Yukawa terms like (jH^ . will be specified 
in the two distinct cases of a particular d = 2 model featuring a column-valued Dirac spinor on the 
codimension, and the generic d > 3 models with SO{d) isospinor Dirac fields, which form square 
matrix arrays. These cases will be presented in the following two subsections, 4.2.1 and 4.2.2. 

In both 4.2.1 and 4.2.2, the essential procedure is to so select the separability Ansatz (jH)). 
such that the resulting residual Dirac equation (fT^ turns out to be a known problem leading to 
normaisable zero modes satisfying (fT^ . 

The analysis of the residual Dirac equation p2|) for these models will be restricted to 'soliton' 
backgrounds of the associated Goldstone models, cf. Appendix C, rather than the corresponding 
backgrounds of the Higgs models, cf. Appendix B. The former are the gauge decoupled versions of 
the latter and the results of these analyses are qualitatively the same. Accordingly, (fT^ is effectively 
replaced in what follows by 

(r'"9^ + /iec)^ = 0. (47) 

i.e. with dm replacing the Dm in (fT^. 

4.2.1 Type (^^) model with /i 7^ : d = 2 

For d = 2 the use of a doublet (0^,0^) with a scalar and a pseudoscalar component yields an 
interesting special situation. We consider f/(l) gauge fields, such that in this model ip{xm) does not 
carry an isotopic index. Specifying the separability Ansatz (jSJ with (jHHj) . namely by 

we end up essentially with the model of jTj ^]. The most interesting feature of this model is the 
presence of fermionic zero modes for Abelian backgrounds with any vorticity n. 

With (jH^ in the residual Dirac equation (fT^. one proceeds to solve the latter in the background 
of the d = 2 Higgs 'soliton', namely the usual Nielsen-Oleson vortex or another {p > 2) member of 
the hierarchy in Appendix B, e.g. the vortex of the system ()B.8|) . 

^When instead the alternative asymptotics pointed out in footnote 2 is adopted, then f2{r) satisfies and it is 
/i(r) that must be rejected. 



13 



Alternatively, as in effect we will, one can solve (|12j) in the background of the d = 2 Goldstone 
'soliton' of the p = 2 member of the hierarchy in Appendix C, namely the vortex of the system ()C.3|) . 
i.e. the one resulting from the gauge decoupling of the system ()B.8|) . We restrict the subsequent 
analysis to that of (fT^ in the associated Goldstone 'soliton' background. 

Substituting (jHHj) with 0™ given by ()28|) and the radially symmetric Ansatz (jHBj) for ip, in (|T^ for 
d = 2, the latter separates for m' = m + 1, and reduces to the pair of coupled first order equations 

rn 

f[--fi + vhf2 = (48a) 
r 

/^ + ^^^^/2 + r//^/i = 0. (48b) 
r 

The Dirac equations in [HlllZlin] reduce to Eqns. ()48aj) . (j48b|l . reproducing the d = 2 result of [310], 
for completeness. 

4.2.2 Type (ii) models with 7^ : d > 2 isospinor ^{Xm) 

The situation here is similar to the case of 'instanton' backgrounds considered in 4.1 and likewise 
the isospinor Dirac field subject to spherical symmetry is 

^ = /i(r)]l + /2(r)r„x„. (49) 

Note here that in ()49p we have in all d dimensions, while in fl45|) we have chiral matrices in 
d = 2n, even, dimensions. 

The separability Ansatz in these cases is specified by (jSH), namely 

e = I ® r^0"^ . 

The residual Dirac equation ()47|1 in the Goldstone 'soliton' background now separates^ and yields 
the following pair of first order equations 



leading to 



f[+Vhfi = 
r 

= Q-V fhdr 



(50a) 
(50b) 
(51a) 



^ = ^e-/^^^ (51b) 
Given the asymptotics of h{r), ()29|1 . and the behaviour of h{r) near the origin to be 

h{r) ^ br , (52) 

^Note here that for the case d = 2, where the Abehan Higgs (or Goldstone) background (|28|l is radial for all vorticity 
n, this separation can take place only for the unit vorticity n — 1 background. This contrasts with the model of |10ll7| 
presented above in 4.2.1. 
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it follws that both fi and /2 vanish asymptotically in the r ^ 1 region as required, but only fi 
converges in the r -C 1 region while /2 diverges and must be rejected. The result is one normalisable 
zero mode, /i (r) . It follows from (j29p that the solution fi (r) has an exponential decay at infinity. 

The corresponding result, for (fT^ in the Higgs 'soliton' background, can be readily found using 
the Ansatz (j221), with the matrices Em'* there replaced by F^, viz. flB.4|) in Appedix B. The result 
is qualitatively the same, with the zero mode fi{r) still localised exponentially, except that the 
energy density of the YMH background brane is power localised rather than that of the Goldstone 
background brane analysed here, which is exponentially localised. 

Instead of one such zero mode, it is possible to construct a family of such solutions by relaxing 
the constraint of spherical symmetry. In particular, imposing only axial symmetry characterised by 
a vortex number n, a family of such solutions labeled by n can be found. We do not present the 
details here. 

5 Summary 

We have addressed the problem of extending the mechanism of confining a fermion to the brane in 
a 4 + 1 dimensional model, proposed in , to the case of 4 + dimensional models, for arbitrary d. 

In the model of ^1 E] , the confinement mechanism relies on the fact that a scalar field model 
in the 1 dimensional extra coordinate, i.e. the codimension-1, supports topologically stable 'soliton' 
solutions. This scalar field enters the 5 dimensional fermionic model through a Yukawa interaction 
term and results in the Dirac equation of the system developing a mass term asymptotically in the 
codimension, which is responsible for the confinement. 

To extend this mechanism to higher dimensions, it seems [2] natural to employ some field the- 
oretic model on the codimension-d > 1 that supports topologically stable finite energy solutions. 
We considered candidates for such models which support either 'soliton' or 'instanton' like solutions. 
Our nomenclature throughout was that 'solitons' are supported by Higgs or Goldstone models, most 
notably featuring dimensionful scalar fields whose self-interaction potential leads to symmetry break- 
ing. 'Instantons' on the other hand are supported by purely YM models or by sigma models, in even 
dimensions. 

We proposed two types of models, in both of which the separation of the Minkowski space 
coordinates from the codimension-c? coordinates Xm was effected by an Ansatz, which also resulted 
in the dimensional descent of the Dirac equation in 4 + c? dimensions, to one in d dimensions, which 
we referred to as the residual Dirac equation. The solutions of the latter were what described the 
localisation of the fermion to the brane. 

The first type, [i), of models was characterised by a Dirac operator, which featured a partial 
derivative in all components of the differential operator. Consequently, the residual Dirac spinors 
were isoscalars. The information on the topologically stable solutions on whose background the 
residual Dirac equation was solved, was encoded in a scalar coefficient in the Yukawa term. This 
quantity was a descendent of the topological invariant of the background system. It was found that 
the residual Dirac equation of these models did not support normalisable zero modes. 

The second type, (ii), of models was characterised by a Dirac operator, which featured a partial 
derivative in the Minkowskian components of the differential operator and a covariant derivative for 
the components on the codimension. Consequently, the residual Dirac spinors were isospinors for 
d > 3 when the gauge group was non Abelian, and only in the d = 2 case when the gauge group 
was Abelian it was isoscalar. Type (n) models did result in normalisable zero modes for the residual 
Dirac equations, provided that the Yukawa term was chosen appropriately, and in the case of pure 
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YM 'instanton' backgrounds this meant its absence. When 'instanton' backgrounds were employed 
locahsation to the brane featured a power of r, while employing 'soliton' backgrounds of Higgs or 
Goldstone models resulted in exponential localisation. 

A nontrivial aspect of our results is that in the models presented in sections 4.2, the Higgs models, 
i.e. the YMH systems, can be contracted down to the associated Goldstone models by elimination of 
the gauge fields, cf. Appendix C. Whether this additional feature in constructing such models is of 
any physical advantage is not obvious, but we note that the energy density of the Higgs model brane 
system is power localised, while that in the Goldstone case is exponentially localised. Also, in some 
of the work of |7j, a Goldstone model has been employed, albeit a model with divergent energy. 
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A Yang— Mills models in d = 2n dimensions 

In all even d = 2n dimensions, it is possible to construct Yang-Mills (YM) models which support 
'instanton' like solutions, which are topologically stable and their energy integrals are finite. The 
fundamental relation that ensures the existence of 'instanton' in this hierarchy ^| of YM systems are 
the inequalities stating the lower bounds on the energy integrals, given by the appropriate topological 
charges, namely the Chern-Pontryagin (C-P) numbers. 

Using the convenient notation for the 2p-form p-fold totally antisymmetrised product of the 
curvature 2-form F = F{2), 

-^(2p)/ii/i2---/i2p ~ ^l^il^l2■^tJ■3^J■3■■■■^^J.2p-l^J.2p] 5 totally antisymmetrised in [/ii/i2---/^2p ], (A.l) 
the simplest such inequality for a 2{p + q) dimensional YM system states 

Tr f|F(2g)|2 + /t2(P-5)-[^|F(2p)|2^ >2/€P-'?TrFAFA... AF, p + g times (A.2) 
V (2p)! / 

where k is a constant with the dimensions of length (if p > q). The left hand side defines the energy 
density of the YM system and the right hand side is proportinal to the [p + g)-th C-P density. 

In the special case where p = q, the 4p dimensional YM systems are scale invariant and the 
inequality ()A.2|1 can be saturated. If the gauge group is chosen to be SO^ip) and the gauge fields 
are in the chiral representations thereof then instanton solutions can be evaluated explicitly in the 
spherically symmetric cases, but these do not interest us here. 

What is relevant in the present context is the family of YM models in 4p dimensions from which 
the d{< 4p) dimensional Higgs models described in Appendix B below are constructed. Also relevant 
are the 2{p+q) = 2n dimensional models defined by ()A.2|) . in any even dimension, whose solutions 
provide the 'instanton' backgrounds used in section 4.1. Solutions to both these types of YM models, 
scale invariant or otherwise, are 'instantons' in the sense that at infinity the gauge connection is pure 
gauge satisfying 

Am — ^ g'^dmg , (A.3) 

r— >oo 

such that in the spherically symmetric case (j22|) the 'instanton' profile (j23j) obtains. 
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B Higgs models in d dimensions 



In this Appendix, we describe the Higgs models ^Tj in d dimensions which support finite 

energy topologically invariant sohton solutions. 

In any given dimension d, a hierarchy of Higgs models supporting solitons can be systematically 
constructed by subjecting the p-th member of the Yang-Mills hierarchy jE], c/ Appendix A, in 
dimension Ap > d, to dimensional reduction down to d dimensions. The descent mechanism essentially 
consists of the imposition of a symmetry, which results in the breaking of the gauge group of the 
original Ap dimensional YM system, and at the same time the components of the gauge connection 
on the extra 4p — dimensions appear as Higgs fields in the residual d dimensional system, namely 
in a Higgs model. By choosing the gauge group of the 4p dimensional system suitably, the gauge 
group of the 4p dimensional YM system breaks down to SO{d), yielding the required d dimensional 
SO{d) Higgs model fTj^J^]. For simplicity we will restrict to the scale invariant YM systems for 
the present purpose. 

Now the action density of the Ap dimensional scale invariant YM system is bounded from below by 
the 2p-th Chern-Pontryagin (C-P) density. It turns out that under the descent mechanism described 
in the previous paragraph, this topological lower bound translates to a new lower bound on the 
energy density of the residual d dimensional Higgs model. The lower bound is given by the residual 
C-P density, which now depends on the residual gauge group, not necessarily SO{d). Such lower 
bounds might be described as BogomoVnyi hounds. That this residual C-P density is a topological 
charge density follows from the fact that it can be shown to be a total divergence |27j. whose resulting 
surface integral turns out to be finite subject to the usual symmetry breaking asymptotics of the Higgs 
field provided that the residual gauge connection also exhibits the requisite asymptotics. 

In the particular case of interest, namely when the residual gauge group is arranged to be SO{d), 
the Higgs multiplet is 

$ = r^ (B.i) 

where the index m = 1,2, ..,d labels also the coordinates Xm, in the same notation as above. The 
symmetry breaking condition of the Higgs field can then be stated as 

, (B.2) 



+00 



where t] is the VEV, related to the compactification scale used in the dimensional descent from Ap 
dimensions. For d > 3, i.e. when the residual gauge group is non Abelian, the Higgs field points 
along the unit vector Xm in the r ^ 1 asymptotic region, i.e. on the d — 1 sphere. In this region the 
gauge group breaks down to SO{d — 1) and the connection field decays as r~^. In the Dirac gauge, 
where the Higgs field in the r ^ 1 asymptotic region points along the d-th direction, the connection 
develops a semi-infinite line singularity in the Xd direction, which is an artefact of this gauge. This 
analogy with the familiar case of the monopole ^3] in the d = 3 case is complete and the residual 
SO{d) connection behaves as 

A^^^lg-'dmg, (B.3) 

namely as half a pure-gauge, rather than as one pure-gauge like an instanton. It is for this reason 
that above, we have called the finite energy topologically stable solutions of Higgs models 'solitons', 
in contrast with the corresponding solutions of even dimensional sigma models and YM systems, as 
'instantons'. 

In terms of the spehrically symmetric Ansatz, which is (j22|) with Tjmn now replaced by Pmn 



17 



the asyptotics of the function w{r) corresponding to ()B.3|) are 



\im w{r) 



+ 1 



hm w(r) = . 



(B.5) 



The only exception to the property ()B.3|) is the d = 2 Higgs model, in which case the boundary of 
the space is not sufficiently large to accommodate a Dirac gauge. 

The most familiar Higgs models which can be construced in this scheme descend from the usual 
SU{2) YM system in 4 dimensions, i.e. from the first {p = 1) member of the YM hierarchy [TT] , 
down to d = 3 and d = 2 respectively. In d = 3 one finds the SO (3) Georgi-Glashow model in 
the Prasad- Sommerfield limit, and in d = 2, the familiar Abelian Higgs model. To illustrate this 
scheme further we have to proceed to p > 2, and for the sake of ease of presentation, we restict to 
the first two nontrivial examples. These are the SO{d) Higgs models arising from {p = 2,d = 3) and 
(p = 2,d = 2). 

The d = 3, 5*0(3) Higgs model pH] is defined by the Lagrangian C which is bounded from below 
by the topological charge density g, 



C = Ti(^{Fi,„Dk]<^y + 6\3i{S,F,,} + [Di<^,D,^f 



(B.6) 



3Q6ijkdkTT 
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(B.7) 



in which $ is given by (jB.ll) . S* = r^^ — and the manifestly total divergence form of g is displayed 
in ()B.7|) . Solutions to this system, were constructed in PE] . 

The d = 2, S0{2) or U{1) Higgs model ^\ is defined by the Lagrangian C which is bounded 
from below by the topological charge density g, 



C = A2[(r7^-b|^)F,, + zD[,^*Z},]^]^ + 24Ai(r/^-|v9|^)lA^r + 18Ao(r/^-|^|^)^ (B.8) 



Q 



r^'A,-3i[r^'-r^'\^\' + -My]^D,v* 



+ k 



(B.9) 

^, and again the topological density 



where we have used the complex valued Higgs field (p = 
g is displayed in manifestly total divergence form. 

It is easy to see that the leading terms making a nonvanishing contribution to the integrals of the 
topological charge densities ()B.7jl and ()B.9|1 . respectively, are the magnetic charge of the monopole 
and the winding number of the Nielsen-Oleson vortex. 



C Goldstone models associated to the Higgs models in d 
dimensions 

In this Appendix we define what we have referred to in the above as the Goldstone models associ- 
ated jini 1211 122] to the Higgs models in d dimensions described in Appendix B above, which also 
support soliton solutions. 
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The aspect of ()B.6|) and ()B.8|) that concerns us here is that both possess gauge decouphng hmits, 
which is a consequence of the fact that when the gauge fields are removed [20] from these densities, the 
remaining density still satisfies the Derrick scaling requirement. Indeed, in this gauge decoupling 
was demonstrated concretely for the numerically constructed solutions: specifically, the Goldstone 
'soliton' in is obtained by gauge decoupling the 'soliton' in 19,. This feature is in stark contrast 
with the same procedure for the usual d = 3 and d = 2 Higgs models arrived at from p = 1 YM, i.e. 
the Georgi-Glashow model and the Abelian Higgs model, in which cases the solitons do not persist 
after gauge decoupling. 

It is thus possible |2Dj to find Goldstone models associated to each d dimensional SO{d) Higgs 
model descended from all p > 2 members of the YM hierarchy. We demonstrate this prescription in 
the two examples considered explicitly in Appendix B. 

In the {p = 2,d = 3) model, eliminating the gauge field in ()B.6|1 . ()B.7|1 we find the Lagrangian 
and topological charge density of the associated {p = 2,d = 3) Goldstone model 

C = Tr(^6A3[9,<l',6',<l']2 + 27A2r7^(77'-<l'')'9,<l'2 + 54Ai(r/'-<I>2)4^ (C.l) 



36£j,fc(9fcTr 



- 2r]^<^di<^dj<^ - ^$2(2<I)9i<I) - di<^<^)dj<^ 



(C.2) 



Note that covariant derivatives in ()B.6|) are replaced by partial derivatives in ()G.1|) . This is the 
rationale behind the corresponding replacement of in (fT^ by dm in ()47|) . 

The system described by HC.l^ is related to that considered in except that in the latter, we 
have selected specific values of the couplings A(i), and also replaced the symmetry breaking potential 
in ()(y.ljl by other postitive definite (and symmetry breaking) potentials, without compromising the 
existence of the solutions. Note that in the asymptotic region, ()G.2|) is equivalent to the winding 
number density. 

In the {p = 2,d = 2) model, eliminating the gauge field in ()B.8|) . ()B.9|) we find the Lagrangian 
and topological charge density of the associated {p = 2,d = 2) Goldstone model 

C = X2\d[.,ip*dj]ip\^ + 2AX,ir]'-\^mdM' + l8Xoiv'~\v\Y (C.3) 



Q = —3ieij dj 



(C.4) 



which is the system investigated in [22], except that in the latter two specific symmetry breaking 
potentials in addition to that in ()G.3|) were employed in the numerical construction of the solutions. 

The above described procedure of constructing the associated Goldstone model for any SO{d) 
Higgs model characterised by any {p,d), can be carried out. 
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